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Abstract 

We reanalyze the WA98 single photon data [H at CERN SPS by incorporating several recent 
developments in the study of prompt and thermal photon production from relativistic heavy ion 
collisions |0]. Isospin and shadowing corrected NLO pQCD, along with an optimized scale 
for factorization, fragmentation and renormalization are considered for prompt photon produc- 
tion. Photons from thermal medium are estimated by considering a boost invariant azimuthally 
anisotropic hydrodynamic expansion of the plasma along with a well tested equation of state and 
initial conditions. A quantitative explanation of the data is obtained by combining k x prompt 
with thermal photons, where k is an overall scaling factor. We show that, elliptic flow of thermal 
photons can play a crucial role to distinguish between the 'with' and 'without' phase transition 
scenarios at SPS energy. 



1. Introduction 

The first observation of single photons by the WA98 Collaboration at CERN SPS is con- 
sidered as a well anticipated turning point in the study of relativistic heavy ion collisions using 
electromagnetic probes fljj] . Earlier observations like the one by the WA80 yfl Collaboration, 
provided only a useful upper limit of this study (for recent developments in the field of direct 
photon production from relativistic heavy ion collisions, see Ref. J3, H). The study of electro- 
magnetic radiations, and in particular photons, as a probe of heavy ion collisions is advantageous 
compared to the study of hadrons, for two main reasons. First of all, photons are emitted from 
each and every stage of the expanding system, whereas hadrons are emitted only from the surface 
of freeze-out after suffering strong interactions. Secondly, photons do not suffer final state inter- 
action (for being electromagnetic in nature, their mean free path is larger than the system size) 
and carry undistorted information from the production point to the detector. The major problem 
in the study of single photons from heavy ion collisions arises from the very small signal to 
background ratio. However, recent developments in the background subtraction methods have 
reduced the size of error bars in the direct to decay ratio for photons considerably. We reanalyze 
the single photon WA98 data by incorporating several new improvements in our understanding 
of prompt photon production from heavy ion collisions and considering the latest developments 
in the field of thermal photon production along with a well defined equation of state and suitable 
initial conditions [2]. 
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Figure 1 : [Left panel] Prompt photons from p+p collisions at ys = 19.4 GeV. Experimental results for E704 [6] for p+p 
collisions and those estimated from p+ n C collisions by E629 [7] and NA3 [8] experiments are given for a comparison. 
[Right panel] Effect of isospin and parton shadowing on the production of prompt photons at tJsnn = 17.3 GeV. 



2. Reanalysis of single photon data at CERN SPS 

2.1. Prompt photons 

The study of prompt photon production in p + p collisions has reached a higher level of 
sophistication and all the available data have now been successfully analyzed with NLO pQCD 
treatment [9] without the inclusion of intrinsic kj for protons. In particular, the suppression of 
single photons at large pj for Au+ Au collisions with respect to the single photons resulting from 
p + p collisions at the same nucleon-nucleon center of mass energy may be largely due to the 
difference in valence quark structure of protons and neutrons ifioll . 

We calculate the prompt photon production from p + p collisions using NLO pQCD treatment 
along with an optimized scale for factorization, fragmentation and renormalization (all equal to 
pr/2) at -\fs — 19.4 GeV (without introducing intrinsic kj) and compare our results with 
various experimental data available at that energy (with proper mass number normalization for 
p + 12 C collisions). This comparison is done as no experimental data are available for p + p 
or n + n collisions at the WA98 center of mass energy ( yjs^N - 17.3 GeV) and the available 
data at the closest center of mass energy is at Vs = 19.4 GeV [see left panel of Fig. OX 
We see that the photons originating from fragmentation are about 30% of those produced from 
Compton+annihilation processes. Our result using NLO pQCD matches well with the NA3 [8] 
data, while it underestimates the E704 [6[] and E629 17J] data, same as reported by earlier studies. 

For prompt photon production from 158A GeV Pb+Pb collisions, isospin and shadowing ifTlll 
corrected NLO pQCD treatment is used with the same scaling factor of pr/2 for factorization, 
fragmentation and renormalization. The effect of isospin and shadowing on photon production 
from heavy ion collisions are investigated by calculating nuclear modification factor (Raa) as 
function of pj and xj ( = Ipjf y[s ) for different beam energies. Results for p+p normalized Raa 
as a function of pj for p+n, n+n, and Pb+Pb collisions are shown in right panel of Fig. [TJ We see 
that the photon production from Pb+Pb collisions is suppressed significantly in the intermediate 
and high pj range, compared to the production from p+p collisions. 

2.2. Thermal photons 

For thermal photons, centrality dependent azimuthally anisotropic boost invariant ideal hy- 
drodynamics is used along with different sets of initial parameters. A well tested equation of state 
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Figure 2: Fit to single photon spectra from Pb(158A GeV)+Pb collisions measured by the WA98 [1] experiment for 
To = 0.2 fm/c [left panel] and 1.0 fm/c [right panel] using scaling factor k = 2.7 and 5.9 respectively. 

is used considering a first order phase transition from the plasma state to the hadronic phase at a 
transition temperature T c (~ 164) MeV. The initial energy density profile is taken as proportional 
to the number of wounded nucleons and five different values of tq are considered ranging from 
0.2 fm/c to 1.0 fm/c (in steps of 0.2 fm/c) keeping the total entropy of the system fixed. The time 
evolution of average energy density (e), average temperature (T) and average radial flow velocity 
(vt/c) at different to are compared. We find that the values of (e) (~ T 4 ) changes significantly at 
large r with changing values of to, whereas (T) and (vt/c) are not affected much [see Ref. J2l for 
detail]. Also, the effective temperature, T e s = T + vr)/(l - Vr), (or the blue shifted tempera- 
ture) is calculated as function of proper time at different to to see the combined effect of cooling 
and expansion (velocity). Thermal photons at different To are calculated considering standard 
rates (QGP photons from Ref. [12], and HM photons from Ref. fl3ll ) of photon production for 
0-10% most central collisions with freeze-out energy density of about 0.075 GeV/fm 3 . 

We find that the prompt photon production is about 17% of the total yield measured by WA98 
and the thermal photon result is almost similar to prompt photon production at to = 0.4 fm/c. 
We also note that the thermal photons from hadronic phase are not affected significantly with 
changing to. A quantitative description of the WA98 experimental data is obtained by using 
the relation 'Thermal + k x Prompt' where, k is adjusted to reproduce the photon production at 
Pt = 2.55 GeV/c. For all to, a normalization at the same pj (= 2.55 GeV/c) provides a good 
description of the data in the entire pj range. We find that the scaling factors k = 2.7, 4.9, 5.4, 
5.7, and 5.9 for prompt photons at to = 0.2, 0.4, 0.6, 0.8, and 1.0 fm/c respectively provide a 
good quantitative agreement with the WA98 data [shown in Fig. [2). We argue that the factor k 
for prompt photons accounts for the Cronin effect, in the case of nucleus-nucleus collisions, as 
well as a pre-equilibrium contribution which must surely be included when to is large. 

2.3. Elliptic flow of photons and hadrons 

In a potentially interesting observation we show that, one additional experimental result, i.e, 
the elliptic flow for thermal photons [14] could actually distinguish between the different values 
of to. The elliptic flow results for different to along with the hadronic matter contribution for 
158A GeV Pb+Pb collisions at CERN SPS are shown in right panel of Fig.fJ] We note that several 
earlier studies have explained the WA98 data considering only the formation a hot hadronic gas 
in the collision and without the formation of QGP phase [ 15]. The estimation of photon flow at 
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Figure 3: V2 for thermal photons for different To, along with contributions from hadronic matter. [Right panel] voipj) for 
primary p mesons from Pb + Pb collisions having b = 7 fm at SPS energy, for different tq. 



SPS can distinguish between the two scenarios of 'with' and 'without' phase transitions as the 
nature of vi would be completely different in the two cases. We also compute the particle spectra 
and V2(pr) for several hadrons and show explicitly that both the spectra and elliptic flow results 
remain unaffected with changing values of To for hadrons. \'2{pt) for p mesons at different to are 
shown in right panel of Fig. [3] 

In conclusion, we present a quantitative explanation of the WA98 single photon data at CERN 
SPS by incorporating several recent developments in the field of prompt and thermal photon pro- 
duction from heavy ion collisions. Thermal photons at different to along with prompt contribu- 
tion enhanced by a '/c' factor describes the data quite well in the entire pj range. We also show 
that thermal photon v>2 can distinguish between the different to and phase transition scenarios at 
SPS energies. 
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